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at least one sensed parameter. The analysis module also includes an output terminal coupled to the processor and connectable to an 
^ output device. The output terminal transmits a third signal representing the predictive MOE. 



WO 02/060662 A2 lUMmillMUHIHll 



GB, GR, IE, IT, LU, MC, NL, PT, SE, TR), OAPI patent 
(BF, BJ, CF, CG, CI, CM, GA, GN, GQ, GW, ML, MR, 
NE, SN, TD, TG). 

Published: 

— without international search report and to be republished 
upon receipt of that report 



For two-letter codes and other abbreviations, refer to the "Guid- 
ance Notes on Codes and Abbreviations" appearing at the begin- 
ning of each regular issue of the PCT Gazette. 



WO 02/060662 



PCT/US02/02690 



SYSTEM FOR AND METHOD OF PERFORMING EVALUATION 
TECHNIQUES ON A LOG OR ROUND TIMBER 



5 RELATED APPLICATIONS 

This application claims the benefit of U.S. Provisional Application No. 
60/265,252, entitled SYSTEM FOR AND METHOD OF PERFORMING 
NONDESTRUCTIVE EVALUATION TECHNIQUES ON A LOG OR ROUND 
TIMBER, filed on January 31, 2001. 

10 BACKGROUND OF THE INVENTION 

The invention relates to a system for and method of performing nondestructive 
evaluation techniques on a log or round timber and, particularly, a system for and 
method of performing nondestructive evaluation techniques on a log or round timber 
for assessing the stiffness and modulus of elasticity of the log or round timber. 

1 5 Many decades of inappropriate management practices, or lack of management 

altogether, have produced large acreages of forest stands that are overstocked with 
small-diameter trees of mixed species. These stands are typically low in value, and 
the harvestable material will not cover the costs of needed management treatments. A 
specific example is the interior west region of the United States, where 39 million 

20 acres of ponderosa pine-type forest have lost ecological integrity due to major 

changes in vegetative structure and composition. These changes have been caused by 
control of fire in an ecosystem where historically there were frequent, low-intensity 
stand maintenance fires. Exclusion of fire has led to the current conditions where 
these stands are now at high risk of attack by insects, disease, and stand destroying 

25 wildfires. Restoration, either mechanical or through prescribed fires, can cost $150 - 
$500 per acre. It is essential to find cost-effective products that can be produced from 
the materials available in these stands so that needed management operations such as 
thinning can be implemented to improve the condition of these stands. Economical 
and value-added uses for these removals can help offset forest management costs, 

30 provide economic opportunities for many small, forest-based communities, and avoid 

-1- 



WO 02/060662 PCT/US02/02690 

future loss caused by catastrophic wildfires. Among the issues of great concern for 
engineering applications of these removals are the variability and predictability of 
their strength and stiffness. 

SUMMARY OF THE INVENTION 

5 A critical need for addressing this situation is the development of 

nondestructive technologies for evaluating the potential quality of stems and logs 
obtained from trees in such ecosystems. Static bending, transverse vibration, and 
longitudinal stress wave techniques are frequently used to assess the modulus of 
elasticity (MOE) of lumber. Excellent correlations between MOE values obtained 
1 0 from these techniques have been shown to exist. Even greater correlations exist when 
using developed models that allow for the prediction of static bending properties. 

Accordingly, in one embodiment, the invention provides a method of 
evaluating a log. The method includes the acts of determining a measured modulus of 
elasticity (MOE) of the log, measuring a property of the log, and calculating a 
15 modified MOE based at least in part on the measured MOE and the measured 
property. 

In another embodiment, the invention provides an analysis module for 
evaluating a log including at least one input terminal connectable to the at least one 
input device. The at least one input terminal is operable to receive at least one signal 

20 representing at least one measured property of the log and at least one determined 
parameter of the log determined in response to an energy being applied to the log. 
The analysis module further includes a processor coupled to the at least one input 
terminal. The processor determines a predictive modulus of elasticity (MOE) of the 
log based at least in part on the at least one measured property and the at least one 

25 sensed parameter. The analysis module also includes an output terminal coupled to 
the processor and connectable to an output device. The output terminal is operable to 
transmit a third signal representing the predictive MOE. 
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In yet another embodiment, the invention provides a system for evaluating a 
log. The system includes an input device operable to acquire at least one property of 
the log and to generate a first signal representing the at least one property. The 
system also includes a sensor attachable to the log. The sensor is operable to sense a 
5 stress wave propagating through the log and to generate a second signal representing 
at least one parameter of the sensed stress wave. The system further includes an 
analysis module coupled to the input device. The analysis module is operable to 
receive the first and second signals, to determine a predictive modulus of elasticity 
(MOE) based at least in part on the first and second signal, and to generate a third 
10 signal representing the modified modulus of elasticity. The system also includes an 
output device operable to receive the third signal. 

In even yet another embodiment, the invention provides a software program 
for evaluating a log. The software program includes at least one software module 
stored in a computer readable medium. The software module is executable to receive 
1 5 at least one measured property of the log including a diameter of the log, receive at 
least one determined parameter of the log determined in response to an energy being 
applied to the log, calculate a predictive modulus of elasticity (MOE) based at least in 
part on the diameter and the determined parameter, and output the determined 
modulus of elasticity. 

20 Other features and advantages of the invention will become apparent by 

consideration of the detailed description and accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a schematic diagram of a system for evaluating a log embodying the 
invention. 

25 Fig. 2 is a flow chart of a method for evaluating a log embodying the 

invention. 

Fig. 3 is a schematic diagram representing a three-point, static-bending 
technique performed on a log. 
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Fig. 4 is a schematic diagram representing a four-point, static-bending 
technique performed on a log. 

Fig. 5 is a schematic diagram representing a transverse-vibration technique 
performed on a log. 

5 Fig. 6 is a schematic diagram representing an analogy for the transverse- 

vibration technique shown in Fig. 5. 

Fig. 7 is a schematic diagram representing a stress-wave measurement 
technique performed on a log. 

Fig. 8 is a graph representing a theoretical response at the end of the log in 
1 0 response to a propagating stress wave. 

Fig. 9 is a schematic representation of an experimental setup for a stress-wave 
measurement performed on a log. 

Fig. 10 is a schematic representation of an experimental setup for a transverse- 
vibration measurement performed on a log. 

1 5 Fig. 1 1 is a front view of an experimental setup for a static-bending 

measurement performed on a log. 

Fig. 12 is a table representing the physical characteristics of red pine and jack 
pine logs. 

Fig. 13 is a table representing the MOE of red pine and jack pine logs. 

20 Fig. 14 is a table representing the results of linear regression analyses of 

various MOEs of red pine and jack pine logs. 

Fig. 15 is a table representing the results of linear regression analyses of 
various MOEs of red pine and jack pine logs. 
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Fig. 16 is a graph representing the relationships of stress-wave MOE to static 
MOE for red pine logs. 

Fig. 17 is a graph representing the relationships of stress-wave MOE to static 
MOE for jack pine logs. 

5 Fig. 1 8 is a graph representing the relationships of stress-wave MOE to 

vibration MOE for red pine logs. 

Fig. 19 is a graph representing the relationships of stress-wave MOE to 
vibration MOE for jack pine logs. 

Fig. 20 is a graph representing the relationships of vibration MOE to static 
1 0 MOE for red pine logs. 

Fig. 21 is a graph representing the relationships of vibration MOE to static 
MOE for jack pine logs. 

Fig. 22 is a graph representing the relationships of predictive stress-wave 
MOE to static MOE for red pine logs. 

1 5 Fig. 23 is a graph representing the relationships of predictive stress-wave 

MOE to static MOE for jack pine logs. 

Fig. 24 is a table representing the results of regression analysis of flexural 
stiffness (EI) of logs determined by different techniques 

Fig. 25 is a graph representing the relationships of stress-wave EI to static EI. 

20 Fig. 26 is a graph representing the relationships of stress-wave EI to vibration 

EL 

Fig. 27 is a graph representing the relationships of vibration EI to static EL 
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DETAILED DESCRIPTION 

Before any embodiments of the invention are explained, it is to be understood 
that the invention is not limited in its application to the details of construction and the 
arrangement of components set forth in the following description or illustrated in the 
5 following drawings. The invention is capable of other embodiments and of being 

practiced or of being carried out in various ways. Also, it is to be understood that the 
phraseology and terminology used herein is for the purpose of description and should 
not be regarded as limiting. The use of "including," "comprising," or "having" and 
variations thereof herein is meant to encompass the items listed thereafter and 
1 0 equivalents thereof as well as additional items. 

A system for evaluating a log or round timber 100 (collectively referred to 
herein as a "log") is schematically shown in Fig. 1. The system generally includes 
one or more input devices 105, an analysis module 110, and one or more output 
devices 115. The analysis module 1 10 includes a housing 117 having one or more 

15 input terminals 120 and one or more output terminals 125. The terminals 120 and 125 
are connectable to the one or more input devices 105 and the one or more output 
devices 115, respectively. In other embodiments, the analysis module 1 10 is directly 
connected to the one or more input devices 105 and the one or more output devices 
115. In yet other embodiments, one or more devices of the one or more input devices 

20 105 and/or one or more devices of the one or more output devices 1 1 5 are 

incorporated within the analysis module 1 10. Even yet other embodiments and 
arrangements of the one or more input devices 105, the analysis module 110 and the 
one or more output devices 115 will become apparent based on the description below. 

As used herein, the term "connection," and variations thereof (e.g., connect, 
25 connected, connecting, etc.), includes direct and indirect connections. The 

connection, unless specified, may be by mechanical, electrical, chemical, and/or 
electro-magnetic means, or any combination of the foregoing (e.g. electro- 
mechanical). 
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In general, the one or more input devices 105 provide input signals to the 
analysis module 1 10. The input signals include input or data. For example, in one 
embodiment, the input signals include data relating to one or more measured physical 
properties of a log. As used herein, the term physical properties" refers to a 
measured characteristic or trait relating to the material log. Example measured 
properties include a diameter of the log, a radius of the log, a length of the log, a 
cross-sectional area of the log, a volume of the log, etc. 

Another possible input signal includes data or measurements relating to a 
sensed parameter when a force or energy (e.g., a stress wave, a vibration, a 
mechanical displacement force, etc.) is applied to the log. For example, the input 
signal may include information or data resulting from a stress wave (discussed further 
below) being applied to the log, information or data resulting from a transverse 
vibration (discussed further below) being applied to the log, or information or data 
resulting from a force (discussed further below) being applied to the log. For a 
specific example and in some embodiments, the one or more input signals include 
data representing an average time between peaks when a stress wave is applied to the 
log. For another specific example and in other embodiments, the one or more input 
signals include data representing a modules of elasticity for a stress wave 
measurement (MOE S w). Using the input signals, the analysis module 1 10 determines 
a MOE, which may be a predictive or modified MOE, and/or a flexural stiffness for 
the log. The analysis module 1 10 then outputs the resulting determinations) to the 
one or more output devices. 

As shown in Fig. 1, the system 100 includes one or more input devices 105 
that provide one or more input signals to the analysis module 1 10. In one 
embodiment of the invention, the one or more input devices 105 include one or more 
operator-controlled input devices for entering data. The one or more operator- 
controlled input devices allow an operator to control the system 100 and/or to provide 
data to the analysis module 125. Example operator-controlled input devices include 
one or more push buttons, one or more trim knobs, a keyboard, a keypad, a touch 
screen, a pointing device (e.g., a mouse, a trackball), or similar devices. In other 
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embodiments, the one or more input devices 105 include one or more acquisition 
devices. Example acquisition devices include accelerometer, an oscilloscope, a 
Metriguard Model 3 12 Bending Proof Tester, a dial gauge, an electronic recording 
device for obtaining a load-deflection signal, and similar devices. Other input devices 
that may be used with the invention include storage devices for providing data (e.g., 
magnetic storage devices, optical storage devices, etc.), or even a network connection 
that connects the analysis module to another digital device. 

As shown in Fig. 1, the system 100 includes an analysis module 1 10. The 
analysis module 1 10 includes a microprocessor 130, and internal memory 135 having 
program storage memory 140 and data storage memory 145. The analysis module 
110 may include other circuitry (e.g., filters, analog-to-digital converters, drivers, 
power circuitry, etc.), which are not shown, to implement the invention. 

For the embodiment shown, the program storage memory 140 includes one or 
more software modules having instructions, and the microprocessor 130 retrieves, 
interprets, and executes the instructions of the one or more software modules to 
control the system 100. This includes receiving the one or more input signals, 
determining the MOE and/or flexural stiffness of a log using the one or more input 
signals, and outputting the calculation^) to the one or more output devices 1 15. 
Different methods of determining the MOE and/or flexural stiffness are discussed 
below. 

In one embodiment of the invention, the analysis module 1 10 is a standard 
personal computer. In other embodiments, the analysis module 1 10 is a hand-held 
device, a personal data assistant, an Internet appliance or similar device. In yet other 
embodiments, the analysis module 110 includes any number of processors or 
controllers constructed with other analog and/or digital logic circuitry having 
integrated and/or discrete circuit elements. 

As shown in Fig. 1, the system 100 includes the one or more output devices 
115. Example output devices are a video display, a printer, an external storage 
device, a second processor connected via a network, and similar devices. The one or 
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more output devices 115 may be incorporated with the one or input devices 105. For 
example, in one embodiment, the input and output devices 105 and 115 include a 
touch screen. For another example, the input and output device 1 05 and 1 15 include 
an acquisition device controlled completely or in part by the analysis module 110. 

The components of the system 100 will be further understood by reviewing the 
operation of the system 100. As schematically shown in Fig. 2, a log is provided for 
analysis (act 150). At act 155, the analysis module 1 10 receives the physical 
properties of the log from the one or more input devices 1 05. Different methods for 
receiving the physical properties include: manually entering data using the operator- 
controlled input device; acquiring the measurements with an acquisition device, and 
providing an acquisition signal or measured data to the analysis module 1 10; and 
providing data via a storage device or network connection. 

At act 160, the analysis module 110 receives the sensed parameters from the 
one or more input devices 105. Different method for providing the sensed parameters 
include: manually entering data using the operator-controlled input device; acquiring 
the parameter with an acquisition device, and providing an acquisition signal or 
measured data to the analysis module 1 10; and providing data via a storage device or 
network connection. 

At act 165, the analysis module 1 10 determines one or more modulus of 
elasticities, and/or flexural stiflhess for the log. In some embodiments, the 
determined MOE is a predictive or modified MOE. Different methods for calculating 
the MOE and the stiflhess are discussed below. At act 170, the analysis module 1 10 
transmits the obtained data, the one or more MOE(s) and/or the flexural stiffness to 
the one or more output devices 115. 

Different methods for determining the MOE and/or stiflhess are now 
discussed. 
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/. DETERMINING A MOE USING A STATIC-BENDING TECHNIQUE 

Measuring the modulus of elasticity of a log using a static-bending technique 
involves utilizing the load-deflection relationship of a simply supported beam with 
different loading patterns. The analysis module 1 10 determines the static MOE 
(MOEs) using equations derived from the fundamental mechanics of materials. For 
example, Fig. 3 schematically shows a standard bending configuration. A load P is 
applied at the mid-span of the log 200. An acquisition device (e.g., a dial gauge or 
electronic recording equipment) senses a load-deflection, and transmits a load- 
deflection signal or measured deflection data to the analysis module 1 10. 
Additionally, a length of the log 200 is measured (e.g., by hand or by a second 
acquisition device) and is provided to analysis module 1 10. The analysis module 1 10 
then determines a static MOE using data taken from the linear elastic region of the 
load-deflection curve. 

Specifically, the analysis module 110 uses the equation: 



where P is the load within the proportional limit, L is the beam span, 8 is the 
deflection at the mid-span within the proportional limit, and I is the beam moment of 
inertia. 

For another example of a static-bending technique, Fig. 4 schematically shows 
a four-point bending configuration. Loads P/2 are equally applied at two points of the 
log 200. Similar to Fig. 3, the analysis module 1 10 receives a load deflection signal 
or measured deflection data, and a length of the log 200. The static MOE for this 
example is calculated by the following equation: 



0) 



MOE s = 



PL 3 



4S& 
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where P is the load, a is the distance from the end support to the nearest load point, L 
is the beam span, 8 is the mid-span deflection, and I is the beam moment of inertia A 
specific embodiment for determining a static MOE is provided below. 

2. DETERMINING A MOE USING A TRANSVERSE-VIBRATION TECHNIQUE 

5 To illustrate a method of measuring the MOE of a log using a transverse- 

vibration technique (MOEt), an analogy can be drawn between the behavior of a 
vibrating log 205 (Fig. 5) and the vibration of a mass M (Fig. 6) that is attached to a 
weightless spring and internal damping force. As shown in Fig. 6, mass M is 
supported from a rigid body by a weightless spring 210 whose stiffness is denoted by 
10 K. Internal friction or damping is represented by the dashpot/X A forcing function 
equaling P 0 sin cot or zero is applied for forced and free vibration, respectively. When 
mass M is set into vibration, its equation of motion can be expressed as follows: 



(3) M^ y^yKx = P 0 sinarr 

Equation (3) can be solved for either K or Z). A solution for K will lead to an 
15 expression for a MOE T of a beam freely supported at two nodal points: 

(4) MOE T =^L 

T 12.65/g 

and for a MOE T of a beam simply supported at its ends: 

f r 2 M? 



(5) MOE T = 



2MIg 



In Equations (4) and (5), the MOE T is a dynamic modulus of elasticity (lb/in 2 
20 (Pa)), where/ is th e resonant frequency (H z), W is the beam weigh t 0b (kg)), L is the 
beam span ( in. (m)), I is the beam moment of inertia (in 4 (m 4 )), and g is the 
acceleration due to gravity (386 in/s 2 (9.8 m/s 2 )). 
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Examples of various systems for performing a transverse vibration evaluation 
are shown in: JAYNE, B.A., Vibrational properties of wood as indices of quality, 
Forest Prod. J. 9(1 1), 1959, pp. 413-416; KAISERLIK et ah, Stress wave attenuation 
as an indicator of lumber strength, Forest Prod. J. 27(6), 1977, pp. 39-43; 
PELLERIN, R.F., A vibrational approach to nondestructive testing of structural 
lumber, Forest Prod. J. 1 4(3), 1965, pp. 93-101; ROSS et al, Transverse vibration 
nondestructive testing using a personal computer, Res. Pap. FPL-RP-502 Madison, 
WI:U.S. Department of Agriculture, Forest Service, Forest Products Laboratory, 
1991; and ROSS et al., Nondestructive testing for assessing wood members in 
structures: A review, Gen. Tech. Rep. FPL- GTR-70 (Rev.), Madison, WI: U.S. 
Department of Agriculture, Forest Service, Forest Products Laboratory, 1994, p. 40.; 
which are all incorporated herein by reference. Additionally, a specific embodiment 
for determining a MOE T is provided below. 

3. DETERMINING A MOE USING A STRESS-WA VE-PROPAGATlON 
TECHNIQUE. 

To illustrate a method of measuring the MOE of a log using a stress-wave- 
propagation technique (MOEsw), consider application of one-dimensional wave 
theory to a homogeneous viscoelastic bar 215 (Fig. 7). After a mass M hits the end of 
the bar 215, a compressive stress wave 220 is generated in the bar 215 and travels 
from left to right at a speed C. As the wave reaches the free end 225, it is reflected as 
a tension wave and begins traveling back down the bar 21 5. Energy is dissipated as 
the wave travels through the bar 215; therefore, although the speed of the wave 
remains constant, movement of particles diminishes with each successive passing of 
the wave. Eventually all particles of the bar 21 5 come to rest. 

Monitoring the movement of a cross section near the end of such a bar in 
response to a propagating stress wave results in waveforms that include a series of 
equally spaced pulses 230 (Fig. 8) whose magnitude decreases exponentially with 
time. The propagation speed C of such a wave can be determined by coupling 
measurements of the time between pulses At and the length of the bar L by 
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The MOEsw can be computed using C and the mass density of the bar p: 
(7) MOE„,=C*p 

Although this equation was derived for an idealized one-dimensional case, it 
5 has been shown to exist for actual three-dimensional members so long as the length of 
the wave is large relative to the lateral dimensions of the member (i.e., log). 

Examples of various systems for performing a stress-wave evaluation are 
shown in: ROSS et al., Technique for nondestructive evaluation of biologically 
degraded wood. Experimental Tech. 18(5), 1994, pp. 29-32; Ross et al., Relationship 

10 between log and lumber modulus of elasticity, Forest Prod. J. 47(2), 1996, pp. 89-92; 
ROSS et al., A stress wave based approach to NDE of logs for assessing potential 
veneer quality, Part 1. Small-diameter ponderosa pine, Forest Prod. J. 49(1 1/12), 
1999, pp. 60-62; SCHAD et al., Stress wave techniques for determining quality of 
dimensional lumber from switch ties, FPL-RN-0265, USDA Forest Serv., Forest 

15 Prod. Lab., Madison, Wis., 1995; WANG et al., Nondestructive methods of 

evaluating quality of wood in preservative-treated piles, Res. Note FPL-RN-0274, 
Madison, WI: U.S. Department of Agriculture, Forest Service, Forest Products 
Laboratory, 2000, p. 9; and WANG, X., Stress wave-based nondestructive evaluation 
(NDE) methods for wood quality of standing trees, Ph.D. Dissertation, Michigan 

20 Technological University, Houghton, MI, 1999; which are all incorporated herein by 
reference. Additionally, a specific embodiment for determining a MOEsw is provided 
below. 

* DETERMINING A MODIFIED MOE USING A MODELS 

In some embodiments, modules are used to modify or predict a MOE. For 
25 example, a regression model may be used to modify a MOE measured by a stress- 
wave technique to predict what the MOE would be for a static bending technique. An 

example mathematical linear regression model equation is: 

-13- 



WO 02/060662 



PCT/US02/02690 



(8) 



y = a+bx 



where y is the modified MOE, a and b are empirical constants, and x is a measured 
MOE (e.g., MOEt, MOEsw, etc.) by nondestructive evaluation. For example and in 
one embodiment, the analysis module 110 receives a MOE S w from the one or more 
input devices 105, and calculates a predictive MOE using the received MOEsw and 
equation (8). In other embodiments, the analysis module 1 10 receives one or more 
physical properties of the log (e.g., the length of a log) and one or more sensed 
parameters (e.g., At for a stress wave evaluation) from the one or more input devices 
105, and calculates a predictive MOE using equations (6), (7) and (8). The empirical 
constants may be previously stored within the analysis module 110, and may be based 
on experimental testing. Specific examples for determining a predictive MOE using a 
linear regression model are provided below. 

Multivariate regression modules can also be used to predict a MOE. An 
example mathematical multivariate regression model is: 



where y is the predicted MOE ; a, b and c are empirical constants; jcy, is a 
nondestructive MOE (e.g., MOE T , MOEsw, etc.); and x 2 is a physical property of the 
log (e.g., a diameter-to-length ration). For example and in one embodiment, the 
analysis module 110 receives a MOEsw from the one or more input devices 105, and 
calculates a predictive MOE using the received MOEsw and equation (9). In other 
embodiments, the analysis module 110 receives one or more physical properties of the 
log (e.g., the diameter of a log) and one or more sensed parameters (e.g., At for a 
stress wave evaluation) from the one or more input devices 105, and calculates the 
predictive MOE using equations (6), (7) and (8). The empirical constants may be 
previously stored within the analysis module 110, and may be based on experimental 
testing. Specific examples for determining a predictive MOE using a multi-variable 
regression model are provided below. 



(9) 



y = axi b x 2 ' 
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While two equations for detennining a predictive MOE was provided, it is 
envisioned that other equations or relationships may be used to predict a MOE. 

J. DETERMINING A FLEXURAL STIFFNESS 

Of the properties and parameters that can be measured nondestructively, e.g., 
5 density, appearance, MOE, and stiffness, etc., stiffness is used most frequently to 
predict the strength of wood materials. Flexural stiffness (EI) is expressed as the 
product of the moment of inertia (I) and modulus of elasticity (MOE) in bending. For 
logs, the moment of inertia is given by 

(I0) '"i* 

10 where d is the average diameter of a log. Upon obtaining a MOE, the flexural 

stiffness of a log can be easily calculated. Of course, the analysis module may receive 
the measured physical properties and the sensed parameters and directly calculate the 
flexural stiffness. 

An example study was performed for comparing various MOEs and stiffness 
1 5 determinations. A discussion of the study is provided below. 

A. MATERIALS AND METHODS 

First, a sample of small-diameter trees were selected from stands and 
harvested to obtain logs. Physical properties (e.g., diameters, moisture contents, and 
densities) of the logs were then measured. This was followed by a sequence of 
20 nondestructive tests using longitudinal stress wave, transverse vibration, and static 

bending techniques to obtain various MOEs and Els of each log. Statistical analyses 
were then used to examine the relationships between log properties determined by 
different techniques. 

A total of 159 small-diameter logs, including 109 jack pine (Pinus banks iana) 
25 and 50 red pine {Pinus resinosa Aft.)> were nondestructively evaluated in this study. 



-15- 



WO 02/060662 



PCT/US02/02690 



These logs came from trees that were grown on the Ottawa National Forest and the 
Lake Superior State Forest in Northern Michigan. 

The jack pine logs used in this study were obtained from an over-age stand of 
jack pine, which is beginning to show signs of deterioration. Ranger District 
5 personnel are able to visually identify four categories of trees in this type of stand: 

live healthy trees (merchantable live), live trees that are showing signs of being under 
stress (suspect), trees that are dead but still containing merchantable material 
(merchantable dead), and dead trees that have deteriorated to the point of having no 
merchantable material (unmerchantable dead). The forest is treating considerable 

10 acreages of these jack pine stands through commercial salvage sales. To be able to 
properly estimate the value of these stands, better information on the value of each of 
the four categories of trees is needed/ Trees of each of these categories were selected 
for this study to address this need. The estimated ages of these jack pine trees ranged 
from 50 to 70 years old. The diameter at breast height (DBH) of sampled trees ranged 

15 from 5.0 to 12.2 inches (127 to 310 mm). 

Red pine logs were obtained from 38 years old research plots that had stocking 
level as the main treatment. The objective of the original study is to examine the 
growth of red pine over time at various stocking levels and correlate volume yield 
with financial yield at the different initial stocking levels. Plots at five levels of 
20 stocking were available 220, 320, 420, 620, and 820 trees per acre. Ten trees were 
harvested from each of the stocking level plots. The DBH of sampled trees ranged 
from 4.70 to 1 1 .50 inches (1 19 to 292 mm). 

After these sampled trees were harvested, a 16-ft- (4.88 -m-) long butt log was 
bucked from each tree on site and then transported to the Forestry Sciences Lab, 
25 USDA Forest Service, North Central Research Station in Houghton, Michigan for 
various nondestructive testing. Upon arrival at the Forestry Sciences Lab, a 2-ft- 
(0.61-m-) long section from each end of the butt log was then cut off and sent to the 
Forest Products Laboratory at Madison, Wisconsin for pulping studies. The 
remaining 12-ft- (3.66-m-) long logs were then used for the purpose of this study. In 
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order to determine moisture content (MC) of sampled trees, 3 cookies were cut from 
butt, middle, and top of each tree respectively. Green weight and oven-dry weight of 
these cookies were then obtained and used to determine tree MC. For each 12-ft- 
(3.66-m-) long log, the green weight and the diameters of both ends were measured to 
5 obtain the green density and the moment of inertia of the log. 

Each log was first evaluated using a longitudinal stress wave technique to 
obtain an estimate of dynamic modulus of elasticity (MOEsw) of the log. Fig. 9 
shows the experimental setup 300 for the stress wave measurements on the logs. An 
accelerometer 305 was attached to one end of the log. A stress wave 310 was 

1 0 introduced to the log through a hammer impact 3 1 5 on the opposite end, and the 
resulting stress wave 310 was recorded using a personal computer 320. A detailed 
description of the instrumentation and analysis procedures can be found in ROSS et 
al., Technique for nondestructive evaluation of biologically degraded wood, 
Experimental Tech, 18(5), 1994, pp. 29-32, which is incorporated herein by reference. 

15 A discussion of its application to large specimens is included in SCHAD et al., 1995 
Stress wave techniques for determining quality of dimensional lumber from switch 
ties, FPL-RN-0265, USDA Forest Serv., Forest Prod. Lab., Madison, Wis., 1995, 
which is incorporated herein by reference. From stress wave measurements, the stress 
wave speed (C) in a log was determined by equation (6). The dynamic modulus of 

20 elasticity (MOEsw) was then calculated using equation (7). 

After stress wave tests, the logs were vibrated using a transverse vibration 
technique. Fig. 10 shows a schematic representation of the experimental setup used 
for the transverse vibration measurement. A digital oscilloscope 325 and an 
accelerometer 330 were used in this test. The log 335 under test was supported at one 

25 end by a knife-edge support 340, and at the opposite end, by a point support 345. The 
accelerometer 330 was located in the middle of the log 335 and glued on the upside 
surface, where the baric was removed or polished to improve the contact between 
accelerometer 330 and the log 335. The log 335 was then set into excitation by 
impacting the middle part of the log using a rubber hammer 350. The free-vibration- 

30 response of the log was observed in the oscilloscope 325. Note that the signal 
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observed was a series of pulses with a gradually decreasing (i.e., decaying) amplitude. 
The vibrational parameter measured was natural frequency. The value for the 
dynamic modulus of elasticity (MOEy) of the logs was calculated with equation (5). 

Static bending tests were then performed on the logs to obtain the flexural 
5 stiffness (EI) and static modulus of elasticity (MOEg) of the logs. Measuring MOE of 
a member by static bending techniques has been widely considered as the foundation 
of lumber grading and NDE of wood and wood-based materials. However, this 
technique is rarely used to evaluate the MOE of logs as a standard method. 
Consequently no standard testing procedure exists for testing small-diameter logs. 
10 Even though, we assume the MOEs of logs to be the real MOE for logs, and used the 
MOEs to evaluate the dynamic modulus of elasticity of logs determined by stress 
wave and transverse vibration techniques. 

A Metriguaid Model 312 Bending Proof Tester 355 (Fig. 1 la) was used to 
conduct static bending tests on all logs. Fig. lib shows the general bending 

15 configuration used. The testing machine was originally designed for proof-loading 
dimensional lumber. In order to test logs, we modified the two end-supports to fit to 
the geometrical shape of small-diameter logs. The modified supports allow testing 
logs with a maximum diameter of 12 inches (305 mm). The span between two 
supports was set as 1 15.5 inches (2.93 m). The distance from loading point to the 

20 nearest support was 38.5 inches (0.98 m), one third of the span. A load was applied to 
the log through two bearing blocks. Deflection was measured in the central region, a 
zone of pure bending without shear deformation. The log under test was first pre- 
loaded to 100 lbs (445 N) and the deflection was set to zero. This procedure was 
mainly used to improve the contact between log, supporters, and bearing blocks and 

25 eliminate the effect of bark on the deflection measurement. The log was then loaded 
to 0.2 inches (5.08 mm) deflection. The practical load value corresponding to this 
deflection was then recorded. The static modulus of elasticity (MOEs) of the log was 
then calculated using equation (2). 
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A RESULTS AND DISCUSSION 

i. Physical Characteristics 

Table 1 (Fig. 12) summarizes the various physical characteristics of the red 
pine and jack pine logs. The average diameters of the butt logs obtained from the 
5 trees ranged from 4.4 to 10.2 inches (1 12 to 259 mm) for red pine and from 4.67 to 

10.99 inches (1 19 to 279 mm) for jack pine. This is a typical diameter range of small- 
diameter timber. For both species, the average moisture content exceeded the fiber 
saturation point (about 30%). However, red pine logs apparently have a much higher 
MC level than jack pine logs. The individual values ranged from 88.2 to 144.6 

10 percent for red pine logs and from 31.2 to 6S.0 percent for jack pine logs. The low 
MC level for jack pine logs was caused by their particular tree source which included 
live, suspect, and dead trees. The suspect and dead trees had already lost a lot of 
moisture by the time they were harvested. Therefore, the moisture content of some 
logs obtained from dead trees was close to or even lower than the fiber saturation 

15 point. 

It was also noted that red pine logs have higher density than jack pine logs. 
The density values for red pine logs ranged from 48.0 to 56.5 pcf (0.77 to 0.90 g/cm 3 ), 
and those for jack pine logs ranged from 28.66 to 53.73 pcf (0.46 to 0.86 g/cm 3 ). The 
lower value and large range of density for jack pine logs was also due to their 
20 particular tree source. 

In appearance, jack pine logs show differences from red pine logs in terms of 
stem shape in cross section and straightness of logs. Red pine logs are mostly round- 
shaped and very straight Whereas some jack pine logs have more irregular shape 
(not round in cross section) and curved stem, which could introduce errors in the 
25 determination of density and moment of inertia of these logs. 

ii. MOEoflogs 

Results obtained from various NDE measurements of both red pine and jack 

pine logs are summarized in Table 2 (Fig. 13). The basic statistics of dynamic MOE 
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(determined by the stress wave and transverse vibration techniques) and static MOE 
(determined by the static bending technique within elastic region) for both species 
were given in this table. 

The static MOE (MOE,) of logs ranged from 0.45 to 1.21 X10 6 lb/in 2 (3.10 to 
8.34 GPa) with a mean value of 0.80 x 10 6 lb/in 2 (5.52 GPa) for red pine and the 
range for jack pine logs was 0.17 to 1.48 xlO 6 lb/in 2 (1.17 to 10.20 GPa) with a mean 
value of 0.81 xlO 6 lb/in 2 (5.58 GPa). It was found that the stress wave technique 
produced a higher estimate of MOE for both species. For red pine logs, the mean 
MOEsw is 1 1.8% and 18.8% greater than its vibrational and static counterpart 
respectively. For jack pine logs, the mean MOE S w is 21.6% and 24.7% greater than 
its vibrational and static counterpart. We believe that the higher value of MOEsw 
could be related with the wave propagation mechanism, dimension of logs, and the 
moisture state of wood in logs. 

The stress waves traveled faster in the outer portion of the wood (i.e., the 
mature wood) when it was propagated through a log in the longitudinal direction. 
This led to a higher stress wave speed on a log and increased the value of the MOEsw, 
which in turn overestimated the MOEsw of the log. It was also found that the 
diameter-to-length ratio could be a critical factor that may affect the stress wave 
behavior in logs. Quantitative analyses of the overestimation in MOEsw, of logs have 
not been reported. 

Compared with the MOEsw of logs, the dynamic MOE of logs determined 
from the transverse vibration technique (MOE v ) is much closer to static MOE of logs. 
The MOEv of red pine logs ranged from 0.58 to 1 .22 x 10 6 lb/in 2 (4.00 to 8.40 GPa) 
and the range for jack pine logs was 0.25 to 1.47 x 10 6 lb/in 2 (1.72 to 10.14 GPa). 
The mean value of MOE v of logs was about 7 percent greater than the mean MOEs of 
logs for both species. 

iii. MOE relationships. 



-20- 



WO 02/060662 



PCT/DS02/02690 



Statistical analysis procedures were used to examine the relationships between 
the various MOE of red pine and jack pine logs. The results obtained from regression 
analyses are presented in Table 3 (Fig. 14) and Table 4 (Fig. 1 5). 

a. Uni viarable regression models 

The correlations among various MOE could be represented by linear 
regression models (y = a + bx). The results of the comparison between three different 
techniques are reported in terms of correlation coefficients that reflect the possible 
reliability of the method for prediction purposes. The square of the correlation 
coefficient expresses the percentage of the total variability explained by the regression 
line. 

In general, the dynamic MOE (MOE sw and MOE v ) of logs was very closely 
correlated with the static MOE (MOEs) for both red pine and jack pine logs. The 
correlation coefficients were found to be 0.87 (MOEsw vs. MOEs) and 0.97 (MOE v 
and MOEs) for red pine logs. Those for jack pine logs were 0.77 (MOE S w vs. MOE s ) 
and 0.92. The linear regression analyses indicated that the developed regression 
models were statistically significant at the 0.01 confidence level. 

Figs. 16, 17, 18 and 19 show the relationships of dynamic MOE, predicted by 
stress wave technique, to dynamic MOE predicted by transverse vibration technique 
and static bending MOE for both species. It was apparent that the red pine logs 
produced a better correlation (r = 0.87 - 0.88) than the jack pine logs (r = 0.76 - 0.77). 
This could be attributed to the geometrical differences between two species. It was 
evident that irregular shape (not round in cross section and curved in stem) of some 
jack pine logs could introduce errors in diameter measurements, thus cause errors in 
the density and the MOE determinations, especially in the MOEsw determinations. 

It was also noted that the plotted data points were more heavily concentrated 
below the 45° line than above, thus indicating that stress wave technique yields higher 
MOE values than its vibrational and static counterpart. As was discussed earlier, the 
higher value of MOEsw could be caused by several factors such as wave propagation 
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mechanism, log size, and moisture state of wood. Of these factors, log size (diameter 
D and length L) seems more important because it could affect stress wave behavior in 
logs. A high diameter-to-length ratio (D/L) could cause significant changes in wave 
propagation path with respect to longitudinal direction when applying longitudinal 
stress wave measurements on logs. Therefore, it seems that the effect of log size 
should not be neglected in the MOE regression models. 

The relationships between MOE v and MOE s of red pine and jack pine logs are 
shown in Figs. 20 and 21 . The results from transverse vibration tests demonstrated 
that a significant improvement had been achieved in comparison to the results from 
stress wave tests. In regard to the relationship of MOE v to MOEs, it was found that 
two species could be combined and represented as a single population. The 
correlation coefficient relating the MOE v to MOEs was 0.93. This value indicates 
clearly that the variation caused by these two species does not affect the relationship. 

b. Multi-variable regression models 

In regard to stress wave MOE (MOE S w), it was found that the univariable 
linear regression models resulted in a correlation coefficient of 0.77 - 0.87 with static 
bending MOE. This value, although significant, indicates a greater scatter of points 
about the regression line than had occurred in MOEv. hi an effort to obtain a better 
prediction model for modulus of elasticity of logs, a multi-variable regression model 
relating the static MOE (MOE S w) to the stress wave MOE (MOEsw) and the 
diameter-to-length ratio was developed. The mathematical regression models used in 
this analysis were assumed to be of the form shown in equation (9). 

The MOE of logs predicted by this equation was then compared against the 
static bending MOEs of logs. Results of the regression analyses and values for the 
constants in the equations are summarized in Table 4 (Fig. 15). 

The relationship between stress-wave-predicted MOE using the developed 
multi-variable model and the static MOE of logs are shown in Figure 22 and 23. 
They indicate that a strong relationship exists between stress-wave-predicted MOE 
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and static MOE. Compared with the univariable linear regression model, a significant 
improvement was achieved in the multi-variable models. The correlation coefficient r 
increased from 0.87 (red pine) and 0.77 (jack pine) for the univariable model to 0.95 
(red pine) and 0.86 (jack pine) for the multi-variable model. This showed that the 
diameter-to-length ratio (D/L) had an interactive effect that contributed significantly 
when used in conjunction with MOEsw* 

iv. Flexural stiffness relationships. 

Of the parameters that can be measured nondestructively, e.g., density, 
appearance, MOE, and stiffness, etc., stiffness is used most frequently to predict the 
strength of wood materials. Therefore, it is important to know the relationships 
between the stiffness determined by these three techniques. 

Flexural stiffness is expressed as the product of the moment of inertia (I) and 
modulus of elasticity (MOE) in bending. For logs, the moment of inertia is given by 
equation (10). Knowing the modulus of elasticity of logs determined by these 
techniques, the various flexural stiffness of logs can be easily calculated. 

The relationships between various log stiffness (stress wave EI, vibration EI, 
and static EI) are shown in Table 5 (Fig. 24) and Figs. 25, 26 and 27. It was 
recognized that red pine and jack pine logs show no distinction in regard to the 
stiffness relationship. Therefore, we combined these two species together and treated 
them as a single population. 

The results revealed that the correlations between these nondestructively 
determined stiffness were extraordinarily strong. In Figs. 25 and 26, the stiffness 
from transverse vibration and static bending tests were plotted as a function of stress- 
wave-determined stiffness. It is important to note that, compared with the MOE 
relationships, the correlations between stress wave technique and the transverse 
vibration and static bending techniques have been improved significantly in term of 
flexural stiffness. Regression analyses indicate that a second-order polynomial 
regression model (y=a+bx+cx 2 ) could best fit the experimental data. The correlation 
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coefficients were 0.97 (stress wave EI vs. vibration El) and 0.94 (stress wave EI vs. 

c 

static EI) respectively. The developed regression models accounted for 97 and 94 
percent of observed behavior. 

Fig. 27 shows the comparison between transverse vibration and static bending 
5 techniques in term of flexural stiffness of logs. Just as in the case of the MOE 

relationship, a linear regression model was found to be the best fitting function to the 
experimental data. The correlation coefficient was 0.94, indicating that 94 percent of 
observed behavior has been accounted for. 

C. CONCLUSION 

10 Based on the results of these experiments, it can be concluded that small- 

diameter red pine and jack pine logs can be successfully evaluated by longitudinal 
stress wave, transverse vibration, or static bending techniques. The dynamic MOE 
(MOEsw and MOE v ) of logs was found to be well-correlated with the static MOE for 
both species. 

1 5 Further, it was found that the a physical parameter (e.g., a diameter-to-length 

ratio (D/L)) had an interactive effect that contributed significantly when used in 
conjunction with measured parameters (e.g., a MOEsw). For example, the developed 
multi-variable model relating static MOE (MOEs) to stress wave MOE (MOEsw) and 
diameter-to-length ratio (D/L) was found to be a better predictor for static modulus of 

20 elasticity of logs. Red pine and jack pine logs therefore could be combined and 
represented as a single population in the prediction model. 

Extraordinarily strong relationships were found between various 
nondestiuctively determined log stiffness. Compared with the MOE relationships, the 
correlations between the stress wave technique and the transverse vibration and static 
25 bending techniques were improved significantly in term of flexural stiffness. 

Thus, the invention provides, among other things, a new and useful system for 
and method of performing nondestructive evaluation techniques on a log or round 
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timber. Various features and advantages of the invention are set forth in the following 
claims. 
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CLAIMS 

What is claimed is: 

1 . A method of evaluating a log comprising the acts of: 
determining a measured modulus of elasticity (MOE) of the log; 

5 measuring a property of the log; and 

calculating a modified MOE based at least in part on the measured MOE and 
the measured property. 

2. A method as set forth in claim 1 wherein the act of determining a measured 
10 MOE includes the act of measuring the MOE of the log by a nondestructive 

evaluation technique. 

3. A method as set forth in claim 1 wherein the act of determining a measured 
MOE includes the act of measuring the MOE of the log by a stress wave propagation 

15 technique. 

4. A method as set forth in claim 3 wherein the act of measuring the MOE by a 
stress wave propagation technique includes the acts of: 

attaching a sensor to the log; 
20 introducing a stress to the log that results in a stress wave having pulses; 

measuring a time between two consecutive pulses of the stress wave; 
calculating the speed of the stress wave based at least in part on the measured 
time; and 

calculating the measured MOE based at least in part on the calculated speed of 
25 the stress wave. 
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5. A method as set forth in claim 1 wherein the act of measuring a property 
includes measuring a diameter of the log (n) times, and 

wherein the act of calculating a modified MOE includes the acts of 
determining an average diameter of the log based on the (n) measurements and 
5 calculating a modified MOE based at least in part on the measured MOE and the 
average diameter. 

6. A method as set forth in claim 1 wherein the act of measuring a property 
includes the acts of measuring a diameter of the log and measuring a length of the log, 

10 wherein the act of calculating a modified MOE includes the act of calculating 

a modified MOE based at least in part on the measured MOE, the measured diameter 
and the measured length. 

7. A method as set forth in claim 1 wherein the act of measuring a property 

1 5 includes the acts of measuring a diameter of the log and measuring a length of the log, 
wherein the act of calculating a modified MOE includes the act of determining 
a diameter-to-length ratio of the log based at least in part on the measured diameter 
and the measured length, and calculating a modified MOE based at least in part on the 
measured MOE and the diameter-to-length ratio. 

20 

8. A method as set forth in claim 7 wherein the act of calculating a modified 
MOE includes the act of calculating a modified MOE using the formula 

y = ax*x c 2 , 

where y is the modified MOE, X] is the measured MOE, x 2 is a diameter-to-length 
25 ratio, and a, b and c are constants. 

9. A method as set forth in claim 1 and further comprising the act of calculating a 
flexural stiffness of the log based at least in part on the modified MOE. 
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10. An analysis module for evaluating a log comprising: 

at least one input terminal connectable to the at least one input device, the at 
least one input terminal being operable to receive at least one signal representing at 
least one measured property of the log and at least one determined parameter of the 
log determined in response to an energy being applied to the log; 

a processor coupled to the at least one input terminal, the processor 
determining a predictive modulus of elasticity (MOE) of the log based at least in part 
on the at least one measured property and the at least one sensed parameter; and 

an output terminal coupled to the processor and connectable to an output 
device, the output terminal being operable to transmit a third signal representing the 
predictive MOE. 

11. An analysis module as set forth in claim 1 0 wherein the at least one input 
device includes an operator-controlled input device. 

12. An analysis module as set forth in claim 10 wherein the at least one input 
device includes a data-acquisition device. 

13. An analysis module as set forth in claim 10 wherein the at least one input 
device includes a storage device. 

14. An analysis module as set forth in claim 10 wherein the at least one measured 
property includes a diameter of the log, and 

wherein the processor determines the predictive MOE based at least in part on 
the diameter of the log. 

15. An analysis module as set forth in claim 10 wherein the at least one measured 
property includes a diameter-to-length ratio of the log, and 

wherein the processor determines the predictive MOE based at least in part on 
the diameter-to-length ratio of the log. 
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16. An analysis module as set forth in claim 10 wherein the at least one 
determined parameter includes a measured modulus of elasticity (MOE), and 

wherein the processor determines the predictive MOE based at least in part on 
the measure MOE. 

5 

1 7. An analysis module as set forth in claim 1 6 wherein the measured MOE 
represents a MOE determined in response to a stress wave being applied to the log. 

18. An analysis module as set forth in claim 10 wherein the at least one 

1 0 determined parameter includes an average time between pulses determined in 
response to a stress wave being applied to the log, 

wherein the processor determines the predictive MOE based at least in part on 
the measure MOE. 

15 19. An analysis module as set forth in claim 1 0 wherein the processor Anther 
determines a flexural stiffness base at least in part on the at least one measured 
property and the at least sensed property. 

20. A system for evaluating a log comprising: 
20 an input, device operable to acquire at least one property of the log and to 

generate a first signal representing the at least one property; 

a sensor attachable to the log, the sensor being operable to sense a stress wave 

propagating through the log and to generate a second signal representing at least one 

parameter of the sensed stress wave; 
25 an analysis module coupled to the input device, the analysis module being 

operable to receive the first and second signals, to determine a predictive modulus of 

elasticity (MOE) based at least in part on the first and second signal, and to generate a 

third signal representing the modified modulus of elasticity; and 
an output device operable to receive the third signal. 

30 
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21. A system as set forth in claim 20 wherein the first input device is an operator- 
controlled input device. 

22. A system as set forth in claim 20 wherein the at least one property includes a 
diameter of the log, and 

wherein the analysis module determines the predictive MOE based at least in 
part on the diameter of the log. 

23. A system as set forth in claim 22 wherein the at least one property includes a 
length of the log, and 

wherein the analysis module determines the predictive MOE based at least in 
part on a diameter-to-length ratio of the log. 

24. A system as set forth in claim 23 wherein the at least one parameter is the 
sensed stress wave, 

wherein the sensed stress wave includes at least two pulses, 

wherein the analysis module is further operable to analyze the sensed stress 

wave to determine an average time between the pulses, and 

wherein the analysis module determines the predictive MOE based at least in 

part on the average time between the pulses. 

25. A system as set forth in claim 23 wherein the sensed stress wave includes at 
least two pulses, 

wherein the at least one parameter is an average time between the pulses, 
wherein the analysis module determines the predictive MOE based at least in 
part on the average time between the pulses. 
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26. A software program for evaluating a log, the software program including at 
least one software module stored in a computer readable medium, the software 
module being executable to: 

receive at least one measured property of the log including a diameter of the 

log, 

receive at least one determined parameter of the log determined in response to 
an energy being applied to the log, 

calculate a predictive modulus of elasticity (MOE) based at least in part on the 
diameter and the determined parameter, and 

output the determined modulus of elasticity. 

27. A software program as set forth in claim 26 wherein the software module is 
executable to calculate the predictive MOE by being further executable to 

calculate a measured MOE based at least in part on the at least one parameter, 

and 

calculate the predictive MOE using the measured MOE. 

28. A software program as set forth claim 27 wherein the software module is 
executable to calculate a predictive MOE by being further executable to calculate a 
predictive MOE using the formula 

y = ax?x c 29 

where y is the modified MOE, x, is the measured MOE, x 2 is a diameter-to-length 
ratio, and a, b and c are constants. 

29. A software program as set forth in claim 26 wherein the software module is 
executable to receive the at least one measure property by being further executable to 
receive a length of the log, and 

wherein the software module is executable to calculate a predictive MOE by 
being further executable to calculate the predictive MOE based at least in part on a 
diameter-to-length ratio of the log. 
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30. A software program as set forth in claim 26 wherein the software module is 
further executable to: 

calculate a flexural stiffness of the log based at least in part on the diameter 
and the determined parameter, and 
5 output the calculated flexural stiffness. 
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